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O ! Abstract 

(N 

We present three new approximation algorithms with improved constant ratios for selecting 
n points in n disks such that the minimum pairwise distance among the points is maximized. 

\^ . (1) A very simple 0(n log n)-time algorithm with ratio 0.511 for disjoint unit disks. 

. (2) An LP-based algorithm with ratio 0.707 for disjoint disks of arbitrary radii that uses a 

linear number of variables and constraints, and runs in polynomial time. 

CiD ' (3) A hybrid algorithm with ratio either 0.4487 or 0.4674 for (not necessarily disjoint) unit 

. disks that uses an algorithm of Cabello in combination with either the simple 0(n log n)- 

fyi2 ' time algorithm or the LP-based algorithm. 

O '. 

The LP algorithm can be extended for disjoint balls of arbitrary radii in M'*, for any (fixed) 
dimension d, while preserving the features of the planar algorithm. The algorithm introduces 
a novel technique which combines linear programming and projections for approximating Eu- 
' clidean distances. The previous best approximation ratio for dispersion in disjoint disks, even 

CNJ . when all disks have the same radius, was 1/2. Our results give a partial answer to an open 

question raised by Cabello, who asked whether the ratio 1/2 could be improved. 

04 ' Keywords: Dispersion problem, linear programming, approximation algorithm. 
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Let 7^ be a family of n subsets of a metric space. The problem of dispersion in TZ is that of 
^ ' selecting n points, one in each subset, such that the minimum inter-point distance is maximized. 

^ This dispersion problem was introduced by Fiala et al. [8j as "systems of distant representatives" , 

generalizing the classic problem "systems of distinct representatives". An especially interesting 
version of the dispersion problem, which has natural applications to wireless networking and map 
labeling, is in a geometric setting where 7?. is a set of disks in the plane. 

Let P be a set of n disks in the plane. Dispersion in disks is the problem of selecting n 
points, one from each disk in T>, such that the minimum pairwise distance of the selected points is 
maximized. Dispersion in disks is a hard problem. Fiala et al. jSj showed that dispersion in unit 
disks is already NP-hard. It is not difficult to modify their construction, which gives a reduction 
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from Planar-3SAT, to show that dispersion in disjoint unit disks is also NP-hard. Moreover, by 
a slackness argument [9l [10], the same construction also implies that the problem is APX-hard; 
i.e, unless P = NP, the problem does not admit any polynomial-time approximation scheme. On 
the positive side, Cabello [3] presented an approximation algorithm Placement with ratio 3/8 
for dispersion in arbitrary disks, and another approximation algorithm Centers with ratio 1/2 
for dispersion in disjoint disks. Cabello also showed that applying both algorithms Placement 
and Centers and taking the better solution yields an approximation algorithm with ratio about 
0.4465 (1/2.2393) for dispersion in unit disks. 

We first introduce some preliminaries. For two points p = {xp, Up) and q = {xq,yq) in the plane, 
let \pq\ denote the Euclidean distance between them: \pq\ = \J [xp — XqY + {Vp — Uq)^- A unit disk 
is a disk of radius one. The distance between two disks is the distance between their centers; e.g., 
the distance between two tangent disks of radii ri and r2 is ri + r2- Throughout this paper, disjoint 
means interior-disjoint. 

We now review the previous best ^-approximation for dispersion in disjoint disks, which is 
achieved by a naive algorithm Centers that simply selects the centers of the given disks as the 
points [3J. Let D = {Oi, . . . , be a set of n disjoint disks of arbitrary radii in the plane. For 
each i, let r^ be the radius of fij. For i ^ j, let 6ij be the distance between fij and fij. Let 5 be 
the minimum pairwise distance of the disks in D, i.e., 6 = miuj^j 5ij. Let OPT denote an optimal 
solution and CEN denote the solution returned by Centers. We clearly have CEN = 6. Since the 
disks are disjoint, 

ri + rj < 5ij, i ^ j. 

It follows that 

OPT < m.m{6ij + ri + rj) <2 min 5ij = 26. (1) 

Consequently, the algorithm Centers achieves an approximation ratio of 

CEN 6 _ 1 
OPT - 2^ ~ 2 

for disjoint disks of arbitrary radii. Cabello asked whether this trivial ^-approximation can be 
improved, even when all disks are not only disjoint but also have the same radius [3l p. 72]. 

We start with a very simple and efficient algorithm that achieves a ratio better than 1/2 for 
dispersion in disjoint unit disks. 

Theorem 1. There is an 0{nlogn)-time approximation algorithm with ratio 0.511 for dispersion 
in n disjoint unit disks. 

Using linear programming, we then obtain the following substantially better approximation for 
dispersion in disjoint disks of arbitrary radii. 

Theorem 2. There is an LP-hased approximation algorithm, with 0{n) variables and constraints, 
and running in polynomial time, that achieves approximation ratio 0.707, for dispersion inn disjoint 
disks of arbitrary radii. Moreover, the algorithm can be extended for disjoint balls of arbitrary radii 
in W^, for any (fixed) dimension d, while preserving the same features. 

We next improve the 0.4465-approximation for dispersion in (not necessarily disjoint) unit disks 
by a hybrid algorithm that uses the algorithm Placement of Cabello in combination with either 
the simple 0(n log n)-time algorithm in Theorem [1] or the LP-based algorithm in Theorem [2j 
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Theorem 3. In combination with an algorithm of Cahello, the simple 0{n\ogn)-time algorithm 
in Theorem [I] yields an 0{n'^)-time algorithm with ratio 0.4487, and the LP-based algorithm in 
Theorem\^ yields a polynomial-time algorithm with ratio 0.4674, for dispersion in n (not necessarily 
disjoint) unit disks. 

It is likely that our method for proving Theorem [21 which uses projections for approximat- 
ing distances, and linear programming for optimization, is also applicable to other optimization 
problems involving distances. 

Related work. The dispersion problem in disks we study here is related to a few other problems 
in computational geometry. We mention several results that are more closely related to ours: 

1. For labeling n points with n disjoint congruent disks, each point on the boundary of a dis- 
tinct disk, such that radius of the disks is maximized, Jiang et al. [TOj presented a Y^+e' 
approximation algorithm, and proved that the problem is NP-hard to approximate with ratio 
more than jg^- 

2. For packing of n axis-parallel congruent squares (congruent disks in the Lqo metric) in the 
same rectilinear polygon such that the side length of the squares is maximized, Baur and 
Fekete [T] presented a |-approximation algorithm, and proved that the problem is NP-hard 
to approximate with ratio more than y|. A |-approximation algorithm for a related problem 
of packing n unit disks in a rectangle without overlapping an existent set of m unit disks in 
the same rectangle, has been obtained by Benkert et al. [2]. 

3. Given n points in the plane, Demaine et al. [1] considered the problem of moving them to an 
independent set in the unit disk graph metric: that is, each point has to move to a position 
such that all pairwise distances are at least 1, and such that the maximum distance a point 
moved is minimized. They presented an approximation algorithm, which achieves a good 
ratio if the points are initially "far from" an independent set. However the approximation 
ratio becomes unbounded for instances that are "very close to" an independent set. Observe 
that in this problem, the optimum may be arbitrarily small, i.e., arbitrarily close to 0. 

4. Dumitrescu and Jiang [6j obtained hardness results and approximation algorithms for two 
related geometric problems involving movement. The first is a constrained variant of the 
fc-center problem, arising from a geometric client-server problem. The second is the problem 
of moving points towards an independent set, discussed previously. 

2 A simple approximation algorithm for disjoint unit disks 

In this section we present a very simple approximation algorithm Al with ratio 0.511 for disjoint 
unit disks, and thereby prove Theorem [TJ 

The idea of the algorithm is as follows. Recall that 5 is the minimum pairwise distance among 
the unit disks. Let a = a{5) be a positive parameter to be specified; in particular, at the threshold 
distance 5 = 2 for disjoint unit disks, we have (t(2) = 2.0883 . . ., which is only slightly larger than 
5. Consider the distance graph of the unit disks for the parameter a, which has a vertex for each 
disk, and an edge between two vertices if and only if the corresponding disks have distance at most 
a. If there is a vertex of degree at least two in the distance graph, that is, if there is a disk close to 
two other disks, then a packing argument shows that the minimum pairwise distance of any three 
points in the three disks must be small. Thus simply placing the points at the disk centers already 
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achieves a good approximation ratio. Otherwise, every vertex in the distance graph has degree at 
most one, and the edges form a matching. In this case, the disks that are close to each other are 
grouped into pairs. The distance between the two points in each pair can be shghtly increased by 
moving them away from the disk centers, at the cost of possibly decreasing the distances between 
points in different pairs. 

Let P be a set of n (not necessarily disjoint) unit disks in the plane. The algorithm Al consists 
of three steps: 

1. Compute the minimum pairwise distance S of the disks in D, and for each disk, find the two 
disks closest to it. 

2. If the distance from some disk to its second closest disk is at most a = (t{S), return the n 
disk centers as the set of points. Otherwise, proceed to the next step. 

3. Place a point at the center of each disk. Then, for each disk, if the distance from the disk 
to its closest disk is at most a, move the point away from the closest disk for a distance of 
(cr — 6) /A, so that the two points in each close pair of disks are moved in opposite directions; 
we will show that 6 < a < 5 + 4, thus the distance {a — 6) /A is between and 1, and each 
point remains in its own disk. Finally, return the set of points. 

Algorithm analysis. The bottleneck for the running time of the algorithm Al is simply the 
computation of the two closest disks from each disk in step 1, which takes 0(n log n) time [3 
p. 306]. The other two steps of the algorithm can clearly be done in 0{n) time. For the proof of 
the approximation ratio, define the following function /(s) for s > 0: 

fis) = ^(l + s)2 + l/2 + V3(l + s)2-3/l. (2) 

The function /(•) is increasing and /(O) = ^/S. The justification for step 2 of the algorithm Al is 
the following packing lemma. Here the disk with center O is close to two other disks with centers 
P and Q, respectively; see Figured) 




(a) (b) 

Figure 1: (a) A linkage of the five segments AP, BQ,CO,OP,OQ for three points A,B,C in three unit 
disks with centers P, Q, O, respectively, (b) The extreme configuration: A, P, O are coUinear, B, Q, O are 
coUinear, \AP\ = \BQ\ = \CO\ = 1, \OP\ = \OQ\ = s, \AC\ = \BC\ ^ \AB\ = t. 
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Lemma 1. Let A,B,C be three points in three unit disks with centers P,Q,0, respectively. Let 
s = max{|OP|, |0Q|} and t = mm{\ AC\, \BC\, \ AB\} . Then t < f{s). 

Proof. Refer to Figure HJa) for a linkage of the five segments AP, BQ, CO, OP, OQ under the 
length constraints max{|OP|, |0Q|} = s and max{|^P|, |CO|} < 1. Then, by a continuous 

movement argument (with the point O fixed), it follows that the minimum pairwise distance t of 
the three points A, B, C is maximized by the extreme configuration in Figure [T]^b). 

Now consider the extreme configuration in Figure mb). Let 9 = \IA0B. In l\AOB, we have 



i = 2(1 + s)sin6l =^ sinO = — -, 

^ ' 2(1 + s) 



In /\AOC, we have 



t2 = i2 + (i + g)2_2(i + g)cos(7r-e) ^ cos0= (1 + ^)^ 1 



2(1 + s) 



Prom sin^ 9 + cos^ ^ = 1, it follows that 

t^ + {t^-{l + sf-lf = A{l + sf. 
Substitute x = {1 + s)^ and y = t^ , and we have 

y + {y - X - Vf = 4x, 

which solves to 

y + + a;^ + 1 — 2xy — 2y -\-2x = Ax 
- (23; + l)y + - 2x + 1 = 



y = X + 1/2 ± y/3x - 3/4 

Both solutions of y are valid, but since we are deriving an upper bound for t, we take the larger 
solution 

y = X + 1/2 + v^3x -3/4. 

Thus, in this extreme case, we have 



^2 = (1 + s)2 + 1/2 + ^3(1 + 5)2-3/4 = f{s). 

It follows that, in general, t < f{s). □ 
Consider the following equation in a: 

5 a + S 



f{a) 2(5 + 2)- 

The next lemma confirms that a exists and lies in the desired range: 
Lemma 2. There is a unique solution a to Moreover, 5 < a < 5 + 4:. 



(3) 
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Proof. We first show that, for any s > 0, s + 1 < /(s) < s + 2: 

f{s) = y(i7^)27i/27y^(rT^F^^ > v(i + s)2 = s + 1, 



f{s) = ^J{l + s)■^ + 1/2 + V3(l + s)2 -3/4 < y(l + s)2 + 1 + V4(l + s)2 = S + 2. 

Since the function /(s) is increasing in s for s > 0, the left-hand side of ^ is decreasing in a. 
On the other hand, the right-hand side of ^ is increasing in a. If o" < 5, then we would have the 
inequality 

5 6 6 _ 5+6 a+5 

J{a) -J{8)^ J+2 ~ 2(5 + 2) - 2(5 + 2)' 

Similarly, if cr > 5 + 4, then we would have the inequality 

5 ^ 5 5 ^ {5 + ^) + 5 ^ a + 5 



f{a) - f{5 + 4) (5 + 4) + 1 2(5 + 2) " 2(5 + 2) ' 

Therefore, there is a unique solution o" to (l3|), and 5 < a < 5 + 4. □ 

We now analyze the approximation ratio of the algorithm Al. Let ALG be the minimum 
pairwise distance of the points returned by the algorithm. Let OPT be the minimum pairwise 
distance of the optimal set of points. Let 

We next prove that the approximation ratio of the algorithm Al is at least c, namely that ALG > 
c • OPT, by considering two cases: 

• If the algorithm returns the n disk centers as the set of points in step 2, then there is a disk 
such that the distances from the disk to its two closest disks are at most a. By Lemma [H we 
have OPT < f{cr). Since ALG = 5, it follows that 

^ > ^. (5) 
OPT - f{a) ^ ' 

• If the algorithm proceeds to step 3, then the distance from each disk to its second closest disk 
is more than a. If two disks have distance at most cr, then they must be the closest disks of 
each other, and the movements of points in step 3 ensure that their two points have distance 
at least (5 + 2(cr — 5)/4 = ((T + (5)/2. On the other hand, if two disks have distance more than cr, 
then after the movements their two points have distance at least a — 2{a — 5)/4: = (cr + 5)/2. 
Thus ALG >{a + 5)/2. Since OPT < 5 + 2, it follows that 

ALG ^ a + 5 

OPT - 2(5 + 2)' ^ ' 

By ([!]), ([5]), and ([6]), the algorithm Al achieves an approximation ratio of c{5) for 5 > 0. It can 
be verified that c{5) is an increasing function of 5 for 5 > 0. Thus, for dispersion in disjoint unit 
disks, the approximation ratio is 

c(5) > c(2) = 0.5110... , for 5 > 2. 

This completes the proof of Theorem [H 



6 



3 An LP-based approximation algorithm for disjoint disks 



Warm-up: one dimension. As a warm-up exercise, we first study the dispersion problem on the 
line and on a closed curve. In these two settings the problem can be solved exactly in polynomial 
time. 

Proposition 1. There exists a polynomial-time exact algorithm based on linear-programming for 
Dispersion in disjoint intervals on the line and on a closed curve. 

Proof. For the line, the input consists on n interior-disjoint intervals [ai,bi], . . . ,[an,bn], where 
(i) Oi < bi, for i = 1, ... , n, and (ii) bi < Oj+i, for i = 1, . . . , n — 1. Computing an optimal solution 
Xi E [aj,6j], i = 1, . . . ,n, amounts to solving the following linear program with the n variables Xj 
and 3n — 1 constraints: 



maximize z (LPl) 

o-i < Xi <bi 1 < i < n 
Xj+i — Xi > z, 1 < i < n — 1 



subject to 



For a closed curve of length L, the input consists on n interior-disjoint intervals [ai,6i],..., 
[a„,6„], where (i) aj < bi, for z = 1, . . . , n, (ii) bi < Oj+i, for z = 1, . . . , n — 1, and (iii) < ai, 
bfi < L. Computing an optimal solution Xi G [oj, i = 1, . . . , n, amounts to solving the following 
linear program with the n variables Xi and 3n constraints: 



maximize z (LP2) 

{o,i < Xi <bi 1 < i < n 

Xj+i — Xi> z, 1 < i < n — 1 
xi + L - Xn> z □ 

Two dimensions. Next we present and analyze the approximation algorithm A2 and thereby 
prove Theorem [2j We first introduce some definitions and notations. Let f^i, . . . , 17^ be n pairwise 
disjoint disks of radii ri, . . . ,r„, and centers oi, . . . ,On. Let < A < A' < 1 be two parameters. 
We set A = 1/2 and A' = 3/4 (with foresight) in order to maximize the approximation ratio. For 
i = 1, . . . , n, let and be two disks of radii A • and A' • rj, respectively, that are concentric 
with Oj. Note that C C Oj. Let aij € [— 7r/2,7r/2) be the direction (or angle) of the line 
determined by and Oj. For a G [— 7r/2, 7r/2), let £q, be any line of direction a. For two vectors 
u = (tti, ^2), and V = {vi,V2), their dot product is {u,v) = uivi + U2V2. The scalar projection of v 
onto u is the length of the orthogonal projection of the vector v onto u, with a minus sign if the 
direction is opposite. It is given by the formula 

proJ^^^ = — — . (7) 
\u\ 

For two points, p and q, let projQ(p, g) denote the length of the projection of the segment pq onto 
a line of direction a, i.e., onto the vector (cos q, sin q). For simplicity, we write projjj(gj, gj) 
instead of proj„^^, (g^, g^). 

The idea of the algorithm is as follows: Suppose we restrict the feasible region of each point 
Pi from the given disk Qj to the smaller concentric disk of radius A • r^. We then show the 
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existence of a good approximation for the dispersion problem constrained to the smaher size disks. 
First, observe that the centers of the original disks 17 j are still in the feasible regions for each of 
the n points. So the ^-approximation that we could easily achieve earlier, is still attainable. For 
instance, setting A = yields the algorithm Centers discussed earlier. Second, observe that if 
A is sufficiently small, then the distance between two points (in two smaller disks) can be well 
approximated by the projection of the segment connecting the two points onto the line connecting 
the centers of the two disks. Enclose each smaller disk in a suitable convex polygon Qj, where 
iJi C Qi C u}[ C ^i- The length of each such projection can be expressed as a linear combination 
of the coordinates of the two points, and we can use linear programming to maximize the smallest 
projection length of an inter-point distance. All the constraints in the dispersion problem will be 
expressed as linear inequalities, at the cost of finding only an approximate solution. We now present 
the technical details. 

We start with a technical lemma that will be used in establishing the approximation ratio of 
Algorithm A2. 

Lemma 3. Let A = 1/2, a G [0, 1], and a G [0, 27r). Then 

1 + aA cos a 1 
— / ^ — 7=- (") 

V 1 + + 2a cos a v 2 

Proof. Since the numerator of the left hand side in ([8]) is non-negative, ([8]) is equivalent to 

1 -|- a^A^ cos^ a + 2a\ cos a > A(l + + 2a cos q). (9) 
After reducing the term 2aAcosQ, (j9]) is equivalent to 

l-Fa^A^cos^a > A + Aa^. (10) 
This follows easily from the following chain of inequalities 

A(l + a2) = < 1 < l + a^A^cos^a, (11) 

as required. □ 

A key fact relating projections to distances is: 

Lemma 4. Let A = 1/2. Consider two disjoint disks Hi and ilj at distance 5ij = \oiOj\. Letpi £ flj 
and pj £ be two points. Let qi G he the point on OiPi at distance X\oiPi\ from Oi. Similarly 
define qj G ojj as the point on ojpj at distance X\ojPj\ from Oj. Then 

\PiPj\ ~ \/2' 

Proof. We can assume w.l.o.g. that Oj = (0,0) and Oj = (1,0), so that 6ij = 1. To represent points, 
we use complex numbers in the proof. Let Re(z) denote the real part of a complex number z G C. 
The point pi is represented by Zi, where Zi G C, with \zi\ < rf, hence qi is represented by Azj. 
Since the disks and ^Ij are disjoint, we have ri + rj < 1. The point pj is represented by 1 -|- Zj, 
where zj £ C, with \zj\ < rj; hence qj is represented by 1 -|- Xzj. Write z = Zj — Zi, and note that 
\z\ < \zi\ + \zj\ < ri + rj < 1. Let z = a(cosa -|- isina) be the complex number representation 
of z, where < a < 1, and a G [0,27r). The segments qiqj and pipj are then represented by 
1 + X{zj — Zi) = 1 -|- Az and 1 + z, respectively. 
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Note that 

projjj(q'i, qj) = Re(l + Xz) = 1 + aA cos a, 

and 

\piPj I = 1 1 + z| = 1 1 + a(cos Q + i sin a) \ = \/l + a"^ + 2a cos a . 
Hence by Lemma [3] we have 

projij(g»,gj) _ 1 + aAcosa ^ J_ 
\PiPj\ Vl + + 2a cos a ~ \/2' 

as required. □ 

We now outhne our LP-based algorithm A2 for disjoint disks. Conveniently select disjoint 
convex polygons Qi, i = 1, . . . ,n, such that C C C ilj, for each i = 1, . . . ,n. For instance 
let Qi be an axis-aligned square of side length concentric with cjj. Since Qi are pairwise disjoint, 
Qi are also pairwise disjoint. Moreover, since Qi and Qj are separated by the perpendicular bisector 
of OiOj, for any qi G Qi and qj S Qj, the dot product {qiqj,OiOj) we use in formulating the LP is 
non-negative. We are lead to the following linear program, LPS, with the constraints expressed 
symbolically at this point. A set {qi, . . . ,qn} of n points is sought, where qi = {xi,yi) £ Qi, 
for i = l,...,n. LPS maximizes the minimum pairwise projection on the line connecting the 
corresponding centers of the disks; that is, for each pair the length of the projection of the 

segment connecting the two points qi and qj, on the line connecting the corresponding disk centers 
Oi and Oj. 



maximize z (LPS) 
qi e Qi, I <i <n 



subject to 



projij(gi,gj) > z, l<i<j< 



n 



Writing the linear constraints. Each symbolic constraint qi G Qi is implemented as four 
linear inequalities, one for each side of Qi. Implement each symbolic constraint projjj((7j, qj) > z as 
follows: Let Oj = (^j, Vi) be coordinates of Oi, for i = 1, . . . , n (part of the input). Consider a pair 
where i < j. Recall that Oij £ [— 7r/2, 7r/2) is the angle of the line determined by Oj and Oj. 
Assuming that ■^i < ^2 < • • • < Cnj we have 

- ■ - o^ 
cosoj,- = 7-, smajo = -p p. (IS) 

\OiOj\ \OiOj\ 

Let Oij = (cosajjjSinajj), so that lajjl = 1. Let s77 = {xj — Xi,yj — yi). According to ([7|), 

P^oiijiqi, qj) = ^"^iJl^i) = . ^--^^ = (3.^. _ 3..) ^os aij + (yj - yi) sin aij. 
\aij\ 

As noted earlier the above expression for the projection is always non-negative. Consequently, 
for each pair {i,j), where i < j, generate the constraint: 

{xj - Xi) cos aij + {yj - yi) sin aij > z, 

where cos aij and sin aij are as in (jlSp . 
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Solving the linear program. 

Lemma 5. For any given e > 0, a (1 — e)- approximation of the solution o/LP3 can be obtained in 
polynomial time. 

Proof. Recall that the linear program LPS finds a point set {qi = {xi,yi),i = 1, . . . , n}, qi G Qi C 
Oj, for which the minimum projection is maximized. The constraints of the linear program LPS 
involve irrational numbers, and hence it cannot be claimed that the original LP is solvable in 
polynomial time. However, it is enough to solve the LP up to some precision. For this, it is enough 
to approximate the numbers involved in the constraints up to some precision, which is polynomial 
in the error of the output [12j. Consequently, we can encode each coefficient into a rational number 
with (l/e)*^^^^ bits. Then, for a constant e, each coefficient has a constant number of bits, and the 
LP algorithm runs in polynomial time; e.g., O(n^) or 0(n^'^) using interior point methods. □ 



Establishing the approximation ratio. 

Lemma 6. For any given e > 0, the approximation algorithm A2 can achieve a ratio at least ^—j^ 
for pairwise disjoint disks. 

Proof. Let OPT = C, denote an optimal solution, given by n points pi, . . . ,pn, where pi G such 
that \piPj\ > Ci for all i / J, and \piPj\ = C for a-t least one pair Let qi G be the point on 

Oipi at distance A|ojpj| from Oj. By Lemma [3] and HJ for all i ^ j, 

projjj(g»,gj) _ l + aX cos a ^ J_ 
\PiPj\ Vl + + 2a cos a ~ \f2 

Since \piPj\ > Ci for all ^ 7^ Ji we immediately conclude that 

pioiij{qi,qj) > (14) 

for all i ^ j. 

Since Ui C Qi, by Lemma [5l the point set found by the linear program satisfies qi G Qi G J^j, 
such that 



P^oiijiqi, qj) > ■ C- (15) 



1 -£ 

Since obviously, \qiqj\ > pi'ojjj(gj, g^), we have found points qi £ Qi C 0,i, such that 



Hence the approximation ratio is at least as claimed. For instance, when setting e = 10""^, we 
get a 0.707 approximation. □ 



Reducing the number of constraints to 0{n). Recall that OPT < 25, as shown in ([T]). 
Observe that the LP solution, z* , is bounded from above as (recall that A' = S/4) 

- 4 - 4 
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where (i, j) are a closest pair of disks. It follows that there is no need to write any constraints for 
pairs of disks at distance larger than 75. Indeed, if now is such a pair, the distance between 
two points, one in Qi and one in Qj, is at least 



4 4 



An easy packing argument shows that the number of pairs of disks at distance at most 75 is only 
0(n): this is the same as the number of pairs of points at distance at most 7 times the minimum 
pairwise distance among n points. 

Extension to any (fixed) dimension d. We briefly sketch the differences. The balls Wj and 
a;- are two smaller balls of radii A • rj and A' • concentric with O,, where A = 1/2 and A' = 3/4. 
Qi is any suitable convex polytope in such that C Qj C C Jl, . The proof of the following 
lemma we need, closely mimics the proof of Lemma HI 

Lemma 7. Let A = 1/2. Consider two disjoint balls ^li and ^Ij at distance 6ij = \oiOj\. Let pi G 
and pj G 0,j be two points. Let qi £ LOi he the point on Oipi at distance \\oiPi\ from Oi. Similarly 
define qj G ujj as the point on ojpj at distance X\ojPj\ from Oj. Then 

pmiij{qi,qj) 1 

Proof. We can assume w.l.o.g. that Oj = (0, . . . ,0) and Oj = (1,0, .. . ,0). To represent points, we 
use vectors in the proof. Let = (1,0, ... ,0). Note that |ei| = 1. The point pi is represented by 
Vi, where It'll < r,; hence is represented by Xvl. Since the two balls Q.i and Q-j are disjoint, we 
have ri + rj < 1. The point pj is represented by ei + UJ, where \vj\ < rj; hence qj is represented 
by eT + Al^J. Write v = Vj — vl, and note that \v\ < |W| + |^| < ri + rj < 1. Let a = \v\; clearly 
< a < 1. Let a G [0, 2tt) be the angle made by v with el. We have 

\PiPj\ = \eT + Vj-Vi\ = \ej; + v\. 
\QiQj\ = |er + Xvj - Awl = |eT + Au|. 



By ([7]), we have 



. / X . , {ei,ei+Xv) 

proJi,(gj, qj) = proje5-(ei + Xv) = — = (d, ei + Xv) 

l^il 



= (ei, ei) + (ei, Xv) = 1 + X{ei,v) = 1 + oA cos a, 

that is, the same expression as in the planar case. Further, by the cosine formula we have 

\piPj\ = |er + u| = -\/|eTp + |?7p + 2|ei||U| cos a = \/l + + 2a cos a, 

matching again the expression from the planar case. 
We conclude in the same way by using Lemma [3) 

projjj(g»,gj) _ 1 + aAcosa ^ J_ 
\piPj\ Vl + + 2a cos a ~ \/2' 

as required. □ 
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The symbolic constraints are implemented analogous to those for the planar case. There exists 
a convex polytope Q C and a function f{d) such that uj C Q C u}' C fl, where Q has f{d) 
facets, and w, uj' and are concentric balls of radii 1/2, 3/4 and 1, respectively (note that already 
for d = 5, a concentric unit hyper-cube is not contained in the unit ball!). The polytope Qi is 
a translate of riQ placed at Oj, so that uji C Qi C uj[ C ^li. Each symbolic constraint qi G Qi 
is implemented as f{d) linear inequalities, one for each facet of Qi. Each symbolic constraint 
P^oiij{qi,qj) > z implements the dot products from the proof of Lemma [71 Since (again) there 
is no need to write any constraints for pairs of balls at distance larger than 76, and the number 
of such pairs is linear in n for fixed d p. 211], the total number of constraints is 0{n). The 
approximation ratio remains the same as for the planar case, namely for any given e > 0, e.g., 

0.707 for £ = 10~^. This completes the proof of Theorem [2l 



4 A hybrid algorithm for unit disks 

In this section we prove Theorem[3l For dispersion in (not necessarily disjoint) unit disks, Cabello [3] 
presented a hybrid algorithm that applies two different algorithms Placement and Centers and 
then returns the better solution. We first briefly review Cabello's analysis for his hybrid algorithm, 
then present an improved hybrid algorithm that uses the algorithm Placement in combination 
with either the simple 0(n log n)-time algorithm in Theorem [1] or the LP-based algorithm in Theo- 
rem [2j In the following, let OPT = 2x and 5 = 2fi. We can assume w.l.o.g. that 5 <2, as otherwise 
the unit disks are disjoint. We also record the obvious inequalities: 

5<OPT<5 + 2<4 ^ fj.<x<l + fi<2. (17) 

The algorithm Placement, which runs in 0{ii?) time, achieves a ratio of 



-\/3 + V^x + V3 + 2x 



X 



2 



ciix) = — , forl<x<2, (18) 

Ax 



and a ratio of at least ^ for < x < 1. 

The algorithm Centers achieves a ratio of 



X — 1 

C2{x) = , for X > 1, (19) 

X 

which is at least | for x > 2. 

Since ci(x) is decreasing in x and C2(x) is increasing in x, the minimum approximation ratio of 
the hybrid algorithm occurs at the intersection of the two curves ci(x) and C2(x) for 1 < x < 2: more 
precisely, ci(x) = C2(x) = —. — i = 0.4465 . . . (1/2.2393 . . .) for x = 1 + = 1.8068 .... 

Recall the approximation ratio c{5) of our algorithm Al as a function of the minimum pairwise 
distance 5 of the disks; see Equation (j4]). Now define 

ai(x) = c(2x - 2), forx>l. (20) 

Then we have the following lemma: 

Lemma 8. For any x > 1, the algorithm Al achieves an approximation ratio at least ai(x) for 
dispersion in unit disks. 



12 



Proof. Recall that the function c{6) is increasing in 6. From (jl7p we deduce that 5 > OPT — 2 = 
2(x— 1). Thus our algorithm Al achieves an approximation ratio of at least c{6) > c(2x— 2) = ai{x) 
for X > 1. □ 

It can be verified that for x > 1, ai(x) is an increasing functions of x. If we replace the 
algorithm Centers by our algorithm Al in the hybrid algorithm, then the two curves ci(x) and 
ai{x) intersect at x = 1.7750. . . and, correspondingly, the approximation ratio of the new hybrid 
algorithm is at least 0.4487 . . . (1/2.2284 . . .). 

We next discuss the hybrid algorithm that runs PLACEMENT and A2. Obviously the n disks 
of radius n < 1 concentric with the n input unit disks are pairwise-disjoint. The hybrid algorithm 
runs Placement on the given unit disks and A2 on the disks of radius fi and then returns the 
better solution. Clearly the solution is valid, and it remains to analyze the approximation ratio. We 

will show that for any given e > 0, it can achieve a ratio at least = n _ ^'j . 0.46749 .... 

1+V9-2V6 ^ ' 

By taking then e = 10~^, we get a 0.4674-approximation. 

We start by relating the optimal solution OPT for the unit disks to the optimal solution OPT^ 
for the smaller disjoint disks. 

Lemma 9. For a problem instance with /i E [0, 1], we have OPT^ > OPT — 2(1 — fi). 

Proof. Consider an optimal solution given by n points pi, . . . ,pn, where pi £ Qi, such that \piPj\ > 
OPT, for all i ^ j, and \piPj\ = OPT for at least one pair {i,j). Let qi G be the point on OiPi at 
distance fJ.\oiPi\ from Oj. Obviously, the set {qi : i = 1, . . . , n} is a valid solution for dispersion in 
the disks of radius ^ concentric with the unit disks f^i, . . . , Moreover, since \piqi\ < 1 — /x, for 
i = 1, . . . , n, by the triangle inequality we have \qiqj\ > OPT — 2(1 — n), for all i 7^ j. Consequently, 
OPT^ > OPT - 2(1 - as claimed. □ 

Now define 

X — 1 ~t~ M 

a2{x,n) = for X G [1,2], and G [0, 1]. (21) 

Observe that 02 (x,^) is an increasing function in both arguments. Then we have the following 
analogous lemma for our algorithm A2. 

Lemma 10. Consider a problem instance with fi £ [0,1] and x G [1, 1 + /i] as in ()17p . For any 
given e > 0, the algorithm A2 can achieve an approximation ratio at least [1 — e) ■ a2(x,/^) for 
dispersion in unit disks. 

Proof. By Lemma [H we have OPT^ > OPT — 2(1 — fi). By Lemma [6l for any given e > 0, the 
algorithm A2 can achieve an approximation ratio at least and therefore find a point set in the 
disks of radius with minimum inter-point distance at least 

^•OPT^>^.(OPT-2(l-^)) = ^.(2x-2(l-M)). 

Equivalently, the approximation ratio obtained is at least 

1-e 2x-2(l-^) x-l + ^ 

V2 2x xV2 □ 
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For a fixed ^ € [0, 1], let xi{fi) be the solution in [l,oo) of the equation ci(x) = a2(x,^), or 

(22) 



4rc xV2 
For solving (j22|) , make the substitution y = x — \, and get 



V^y + V4-y2 = 2V2{y + /i) 
(3 - V6)y^ + (4 - ^/6)^2/ + (2;u2 - 1) = 
_ -(4 - ^/6)/i ifc ^Jl2 - 4x/6 - 2//^ 
" 2(3 - ^/6) ■ 

Since the solution with a minus sign is negative, hence infeasible, we get 



I . -(4-x/6)/i+ Vl2-4x/6-2/i2 
yi(/^) = _ > xi(/i) = l + yi(/i). (23) 

Observe that yi(/i) and xi{jx) are both decreasing functions of for < < 1, with yi(0) = 
= 1.34..., and 2/l(l/^/2) = 0. 

V3-V6 

Let /io be the solution in [0, 1] of the equation yi(/i) = /i, or 



2(3 - 

After rearranging terms and squaring, we get 



(4- V6)/x+ Vl2-4^ /6-2/.2 ^24) 



Lemma 11. 

a/2 

ci(l + ^o) = a2(l + /Uo,Mo) = ; ^ = 0.4674.... (26) 

1 + V9 - 2^/6 

Proof. By definitions we have 



ci(l + ;Uo) = ci(l + yi(^o)) = 02(1 + yi(Mo),/^o) = 02(1 + ^0,^0) 
(1 + ^0) - 1 + /^o \/2/io \/2 \/2 



0.4674 . 



^2(^0 + 1) A^o + l 1 + ^ 1 + \/9 - 2V6 

as required. □ 

It suffices to prove the following lemma and then apply Lemma [10] to complete the proof of the 
theorem. 

Lemma 12. For every fi G [0, 1] we have 

min max{ci(x), a2(x, /i)} > : = 0.4674 . . . . (27) 

^e[i,i+M] 1 + y/g _ 2^6 
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Proof. Let h{x,n) = max{ci{x),a2{x, fi)}. We distinguish two cases. 

Case 1: ^ < fio- By the definition of /io we have /io = UiifJ-o) < Observe that 1 < x < 

1 + fi < 1 + /.iQ. By ([26]) . and the monotonicity of ci(-), we have 

\/2 

h{x,n) > ci(x) > ci(l + ^) > ci(l + ^o) = , = 0.4674... . 

1 + V9 - 2^/6 

Case 2: > ^q. If x < 1 + /xq then l<j;<l+/io<l + /U. By (f26l) and the monotonicity of 
ci(-), 

a/2 

^) > ci(x) > ci(l + /io) = , = 0.4674 .... 

1 + V9 - 2^/6 

If X > 1 + /^o then l<l + //o<x<l + /x. By (j26|) and the monotonicity of a2(-, //) and 02 (x, •), 

\/2 

/i(x, ^) > a2(x, /u) > 02(1 + fiQ, > ^2(1 + ^0, Aio) = , ^ = 0.4674 

1 + V9 - 2^/6 



Hence the inequahty (l27|l is satisfied in both cases. □ 
This completes the proof of Theorem [3l 

5 Summary 

Recah our two algorithms Al and A2 and the two algorithms Placement and Centers by 
Cabello [3]. In conclusion, we summarize the current best approximation ratios for the three 
variants of dispersion in disks: 

• Arbitrary (not necessarily unit or disjoint): 3/8 = 0.375 by Placement. 

• Unit (not necessarily disjoint): 0.4674 by A2 plus PLACEMENT, which improves 0.4465 by 
Centers plus Placement. 

• Disjoint (not necessarily unit): 0.707 by A2, which improves 0.5 by Centers. 

Acknowledgment. The authors would like to thank Andreas Bjorklund for pertinent comments 
on an earlier version of the manuscript. 
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